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ABSTRACT 

Using galaxy sample observed by the BATC large-field multi-color sky survey and galaxy data of SDSS 
in the overlapped fields, we study the dependence of the restframe r-band galaxy luminosity function 
on redshift and on large-scale environment. The large-scale environment is defined by isodensity 
contour with density contrast 5p/p. The data set is a composite sample of 69,671 galaxies with 
redshifts 0.03 < z < 0.5 and r < 21.5 mag. The redshifts are composed by three parts: 1) spectroscopic 
rcdshifts in SDSS for local and most luminous galaxies; 2) 20-color photometric redshifts derived from 
BATC and SDSS; 3) 5-color photometric redshifts in SDSS. We find that the faint-end slope a steepens 
slightly from —1.21 at z ~ 0.06 to —1.35 at z ~ 0.4, which is the natural consequence of the hierarchical 
formation of galaxies. The luminosity function also differs with different environments. The value of 
a changes from —1.21 at underdense regions to —1.37 at overdense regions and the corresponding M* 
brightens from —22.26 to —22.64. This suggests that the fraction of faint galaxies is larger in high 
density regions than in low density regions. 

Subject headings: galaxies: distances and redshifts - galaxies: luminosity function - galaxies: evolution 
- cosmology: large-scale environment - cosmology: observations 



1. INTRODUCTION 

Galaxy luminosity function (LF) is a powerful tool in 
the study of galaxy formation and evolution. Galaxy 
LF is directly related to galaxy mass function. Press 
and Schechter (1974) present a simple analytical for- 
mula for the mass distri bution bas e d on the hierarchi- 
cal assembly of galaxies. iSchechteil l)1976l) gives an em- 
perical functional form for galaxy LF. It has been pro- 
posed that it should b e universal l)Luggerlll986t iCollessI 
T989t iTrenthaml fl998^ . However, mass function, star- 
formation process as well as morphological characteris- 
tics of galaxies are affected by their environments and 
evolve with time, galaxy LF is expected to change with 
time and to vary with galaxy characteristics and the den- 
sity environments. Deep wide-area sky surveys, such 
as th e Two-Degree Field Galaxy Redshift Su rvey (2dF- 
GRS MColless et al.ll200lt iNorberg et al.|[200l and Sloan 
Digital Sky Survey CSDSS: IYork^Hdll200TT) . can gen- 
erate large samples with a range of redshifts, which are 
best suited for the mesurment of LF of galaxies. Many 
studies have been done in galaxy LF by morphological 
types HWolf et al. 1 120031 ICross et aljEooi ICroton et al I 



200l. bv redsh i fts lILillv et al J 1199.1 lEllis et, all 119961 



Wolf et, alJl2003l ILovedavl 1200 
vironmcnts llMercurio et all 12 003: 



an d by large-scale en- 
- Hai nes et all 120041: 
2003J). These enable us to 



ICroton et al.l2004HHovle et all 

test the theories of formation and evolvement of galaxies 
i n different cosmolo gical models. 

Ellis et all (Jljj96) indicate that there is indeed a steep- 
ening in the faint-end slope with redshift. Limited by the 
size of data sets, however, the evolution of LF to high red- 
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shift was not constrained well. For galaxy LF by large- 
scale environ ment, analyses reveal that, at low-density 
environment jE fstathio uet alll988l lLovedav et al.ll992t 
IHovle et all 120031. the faint-end slope turns out to be 
a ~ — 1, and at high-density regions the slope seems 
to be steeper with —1-8 < a < —1.3 jDe Pronris et alJ 
119951: iLumsden et al.lll997t iValotto et al.lll997fT 

In this paper, we use our BATC (Beijing- Arizona- 
Taiwan-Connecticut) 15 intermediate-band color sky sur- 
vey data and SDSS 5-color sky survey data to study the 
evolution of LF and the environmental effects on the LF. 
BATC photometric system has an average depth of 20.5 
mag and the corresponding redshifts of most galaxies are 
less than 0.3. The SDSS broadband photometric system 
has an average depth of 23.0 mag w ith most re dshifts less 
than 0.5. From the investigation in lXia et alJ l|2002tl . the 
accuracy of redshift determination by 15 intermediate- 
bands, <j z ~ 0.02, is much better than that of broad- 
band, er 2 ~ 0.05, in the same photometric magnitude 
errors. Therefore, with the accurate photometric red- 
shifts at z < 0.3 and with small — at 0.3 < z < 0.5, 

z 5 

it is possible to combine data of these two systems for 
accurate measurements of galaxy LF to high redshifts. 

The content of this paper is as follows. In § 2 we 
briefly describe the data sample, the application of pho- 
tometric redshift code hyper z and fc-correction. The fit- 
ting of galaxy LF and results of LF evolution are given 
in § 3. Method of environmental classification and re- 
sults of environmental effects on LF are presented in § 
4. § 5 discusses the effects of photometric redshift un- 
certainty on galaxy LF shape by simulation and sum- 
marizes our conclusions. Through out this paper, we 
assume a ACDM cosmological model with matter den- 
sity Q m = 0.3, vacuum density Qa — 0.7, and Hubble 
constant H = 100/ikms _1 Mpc _1 with h = 0.75 for the 
calculation of distances and volumes (see lHogd ll999) . 

2. DATA 
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2.1. Sample 

We use galaxy data including BATC 15-color pho- 
tometries, SDSS spectroscopies and 5-colors photome- 
tries. The BATC Sky Survey performs photometric ob- 
servations with a large field multi-color system. The ob- 
servation is carried out with the 60/90 cm f/3 Schmidt 
Telescope of National Astronomical Observatories, Chi- 
nese Academy of Sciences, (NAOC) located at the Xin- 
glong station. For detailed description of survey and 
performance, see Zhou et al. 2002. The SDSS per- 
forms imaging and spectroscopic surveys over n stera- 
dians in the northern Galactic cap with a 2.5 m tele- 
scop e at Apache Poin t Observatory, Sunspot, New Mex- 
ico ijYork et al.N2000|) . The detailed description of the 
p hotometric and spectros copic parameters can be found 
in lStouehton et alJ l)2002|) . 

We select 17 fields, totally ~ 17 deg 2 , observed by 
BATC and overlapped with SDSS sky survey. Galax- 
ies with the spectroscopic and photometric informa- 
tion in SDSS are obtained from SDSS Data Release 2 
(http:/ /www.sdss.org/dr2/). The data of 69,671 galaxies 
are achieved in SDSS with r < 21.5 and 0.03 < z < 0.5. 
Galaxies in BATC are selected by coordinates given by 
BATC and SDSS. 10,681 galaxies in BATC are obtained 
with the distance deviations in BATC and in SDSS less 
than 2'(0. To combine photometries in BATC and SDSS, 
we need to apply aperture correction to SDSS model 
magnitudes since that an aperture of 4 pixels (i.e., r ap = 
6"8) is ap plied in BATC photometries (see details from 
lYuan et alJl2003|) . The formula of aperture correction is 
as below: 
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J °° 2irrl{r)dr 



(1) 



where m ap is the aperture magnitude, I(r) is the pro- 
file function of surface intensity for the best fit of de 
Vaucouleurs or exponential model. For these common 
galaxies, we estimate photometric redshifts by the total 
20 colors. 

2.2. Redshift 

The redshifts in our catalog are measured by three 
methods. 1,362 galaxies have spectroscopic redshifts ob- 
served by SDSS; 10,681 galaxies have photometric red- 
shifts estimated by 20 color photometries consisting of 
BATC and SDSS; and the rest 57,628 galaxies have pho- 
tometric redshifts estimated by 5 color photometries of 
SDSS. The photometric redshift technique is based on 
SED fitting to estimate redshifts by comparing the spec- 
trum of an object, which should include several strong 
spectral features such as 4000A break, Lyman-forest 
decrement, etc, with the tem plate spectra. We use the 

" to 



decrement, etc, witn tne template spect ra, We use 
hyper z program developed bv lBolzonella et all l)200' 



estimate redshifts. The accuracy of redshift determi- 
nation in BATC has been achieved to be 0.02 and the 
estimated accur acy by broadband filters is about 0.05 
ijXia et alJl2"002]) . 

The accuracy of photometric redshifts is assessed by 
galaxies with spectroscopic redshifts in the fields for both 
systems. Fig. 1 shows the comparision between photo- 
metric redshifts and spectroscopic redshifts. We can dis- 
tinctly sec that redshifts of certain amount galaxies are 
overestimated by the 5-color SDSS photometries. The 



combined 20 colors can take them back to reasonable 
estimation. The redshifts in the range from 0.3 to 0.4 
seem totally to be overestimated in SDSS, and this may 
effect the measure of luminosity function dramatically 
because of the dominant amount of galaxies of SDSS in 
high redshift layer. The uncertainties are a= 0.017 for 
20-color galaxies and a= 0.022 for 5-color galaxies with 
excluding those Az > 0.05. Considered that galaxies 
used to estimate accuracy are bright in magnitudes and 
accurate in photometries, for faint galaxies, the errors 
of estimated photometric redshifts should be larger than 
these estimates. To assess how the uncertainty of photo- 
metric redshift affect the predicted luminosity function, 
we will investigate by simulation in § 5s. Galaxies with 
z < 0.03 are excluded in the construction of galaxy lumi- 
nosity for the large relative errors for local galaxies. Fig. 
2 is the distribution of galaxy redshifts in our compos- 
ite sample. The first histogram with simple line is the 
redshift distribution for the total galaxy sample. The sec- 
ond histogram is that for galaxies with 20 colors. And 
the third histogram with filled area is for galaxies with 
spectroscopic redshifts. Fig. 2 shows us directly the con- 
tributions from these three sources. 
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Fig. 1. — The comparison between photometric and spectroscopic 
redshifts estimated by 20 colors of BATC and SDSS (left panel) 
and 5 colors of SDSS (right panel), respectively. 




Fig. 2. — Redshift distributions for total galaxies (histogram 
with simple line), galaxies with 20-color photometric redshifts and 
galaxies with spectroscopic redshifts (histogram with filled area). 



2.3. k-correction 



Galaxy Luminosity Function at 0.03 < z < 0.5 
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By photometric redshifts we calculate absolute magni- 
tudes in the restframe r-band of SDSS: 

M r = r-25-5logd L (z)-A r -k(z). (2) 
where (Il(z) is the luminosity distance in Mpc, A r is 
the reddening extinction correction due to intergalactic 
and interstellar dust scattering and absorption, and k(z) 
is the fc-correction due to the shift of spectrum by red- 
shift, fc-correction is significant here since that the sam- 
pled redshifts span to high redshifts. The common ap- 
plied method is to estimate fc-correction by galaxy color 
and morphological type. In the advantage of photomet- 
ric redshift fitting, we derive fc-corrections directly from 
the spectra templates. In photometric redshift fitting, a 
best fit template can be achieved for every galaxy, which 
means the most semblable in the spectrum shape. We 
modify slightly the code and output the best fit template 
spectrum at z— 0. Assuming that the observed spectrum 
has the same fc-correction with the best fit template, the 
fc-correction can be corrected directly. The accuracy of 
this correction is investigated by simulation here. 

By the pro cedure make-catalog (see details in 
IXia et alJl200l . we b uild a catalog of 100 galaxies 
with random types of IBruzual et al.l (^993) and red- 
shifts in the range of z < 0.5. Totally 20 filters are 
used. The model fc-corrections are output directly from 
makejzatalog and the estimated corrections are output 
by hyper z. We, hence, can estimate the accuracy of 
the fc-correction. Fig. 3 shows us the distribution of 
the corrections with redshifts. The corrections k(z) in 
the restframe r-band range from to 2 mag. Fig. 4 
shows the comparison between model and estimated cor- 
rections, and the distribution of the deviations. The rms 
error of fc-correction is about 0.05 mag with a small offset 
of 0.01 mag. 

Redde ning extinction is obtained by step fitti ng in 
hyper z i|Xia et, alJ l2002t) . The reddening law of lAllenl 
(1976) for the Milky Way is adopted. The value of red- 
dening correction is that best fitted by the template spec- 
trum. Here, we also assess the accuracy of the total cor- 
rection of reddening extinction A r and fc-correction by 
simulation. It is found that the rms error of this total 
correction is about 0.08 mag with a small offset of 0.01 
mag. In § 5 we will demonstrate by simulation that mag- 
nitude error plays the most important role in the predic- 
tion of luminosity function and this magnitude error will 
effect LF fit slightly. 



Fig. 3. — The distribution of simulated fc-correction with red- 
shifts. 
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Fig. 4. — Assessment of the accuracy of fc-correction and red- 
dening extinction. The left panel is the comparison between model 
correction and estimated correction. The right panel is the distri- 
bution of the deviations. 



3. EVOLUTION OF LUMINOSITY FUNCTION 

Luminosity function is the number density of galaxies 
as a function of luminosity. To measure luminosity func- 
tion, we adopt t he approximation of Schechter function 
<|Schechterlll976l) : 

$(L)dt = 0*(^-rexp(-A)d(il), ( 3 ) 

expressed in the form of per unit absolut magnitude: 

$(M)rfM = 0.41n(10)<M0-°' 4(M - M - )(Q+1) exp[-10-°- 4 ( M - M *)]dM, 

where M* is the characteristic magnitude, i.e., the point 
at which the bright end cutoff sets in, corresponding to 
L*, a is the power-law slope of the faint end, and is 
the normalization constant. We search these parameters 
by x 2 minimization fitting. The fits are performed over 
the magnitude range —24 < Mi < —16 and the normal- 
izations are done in this magnitude range. 

To correct the incompleteness arising from the selec- 
tion effects of distanc e, the traditiona l \/V(M) method 
originally proposed bv lSchmidtJ |l968) is implemented in 
this paper. In the assumption that galaxies distribute 
homogeneously in comoving space, the correction func- 
tion is 1/V(M), where V(M) is the maximum volume 
that determined by the maximum distance at which a 
galaxy with absolute magnitude M can be observed in 
the apparent magnitude limit. The comoving volume 
and luminosity distance are calculated as that given by 
Hogg 1999: 



dVc = 



d L {zf 



d L {z) = (1 + 2)- 



H (l + z) 2 E(z) 
2 dz' 



dfldz, 



(5) 



H J E(z') 



The structure in small-scale and the possible evolution 
in number density with redshfit can produce spurious es- 
timate of local LF in this estimator. For the large size 
of the data set it can partly make up the effect caused 
by inhomogeneity. For the Malmquist Bias in distance 
estimates, it origins from that observation effect that 
greater numbers of galaxies in the univers at greater dis- 
tances and hence more will have been scattered down 
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from larger distances than up from smaller ones. We 
follow the method given bv iLvnden-Bell et alJ l|1988fl to 
correct the distance estimates. The correction formula is 
as below: 

i? = i? e exp[(a+i)A 2 ], (7) 

where R e is the estimated luminosity distance Dl, a — 3 
for uniform distribution, and A is the dispersion in lni? e . 

To study the evolution of galaxy LF, we split our 
galaxy sample into three redshift layers, 0.03 < z < 0.1, 
0.1 < z < 0.3, and 0.3 < z < 0.5, with 5,289, 26,162, 
and 38,220 galaxies respectively. Fig. 5 shows the LFs 
derived for the three layers 0.03 < z < 0.1 (filled trian- 
gles, solid line), 0.1 < z < 0.3 (open circles, dot-dashed 
line), and 0.3 < z < 0.5 (filled circles, dotted line). The 
error bars are the errors of poisson counts. Measured pa- 
rameters are given in Table 1 and the error contours are 
plotted in the left panel of Fig. 6. As shown in Fig. 5 
and Table 1 , the faint-end slope a steepens slightly from 
-1.21±0.02 to -1.25 ±0.03 and -1.35±0.08 with the in- 
crease of redshift. The points appear some discrepancy 
with Schechter function in bright end. The first point 
exceeds than expectation. This may result from several 
reasons. First, it may be due to the observation effects. 
For the nearest galaxies, the uncertain proper motions 
can lead to large errors in rcdshifts and hence overesti- 
mate the most luminous galaxies. On the other hand, the 
amount of intermediate luminous galaxies is much more 
than that of the most luminous galaxies in the obser- 
vation, the uncertainty of magnitude and photometric 
redshift, therefore, can contaminate more intermediate 
redshift galaxies to low redshifts and then bring to the 
excess of most luminous galaxies. 
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Fig. 5. — The luminosity functions fitted for the three redshift 
layers 0.03 < z < 0.1 (filled triangles, solid line), 0.1 < z < 0.3 
(open circles, dot-dashed line) and 0.3 < z < 0.5 (filled circles, 
dotted line). 

For local gal axies 0.03 < z < 0.1, a = -1.21 ± 0.02 is 
consistent with IBlanton et alJ l|200lP a = -1.20 ± 0.03, 
which is achieved by 11 275 galaxies complet e to r < 17.6 
over 140 deg 2 in SPSS. IBlanton et all l)2003ft reevaluated 
galaxy luminosity function at z = 0.1 in the 01 r frame 
by a larger sample of 147,986 galaxies. A much flatter 
faint-end slope a = —1.05 ± 0.01 was found. This large 
difference is due to accounting for the evolution of lumi- 



Table 1 

Shape parameters of LFs fitted to the three 
redshift-binned subsamples and that fitted to high 
and low density subsamples by density contrast 

Sp/p = 30. 



0«(xlO~ 2 ) M, a, x 2 

0.03 < z < 0.1/5289 10.16 ±0.08 -21.80 ±0.16 -1.21 ±0.02 0.43 

O.K z < 0.3/26162 8.56 ±0.13 -21.91 ± 0.12 -1.25 ± 0.03 2.59 

0.3 < z < 0.5/38220 4.42 ±0.17 -22.69 ± 0.21 -1.35 ± 0.08 29.77 

high/54267 4.14 ±0.09 -22.64 ±0.18 -1.37 ±0.04 10.78 

low/15404 8.83 ±0.05 -22.26 ± 0.16 -1.21 ± 0.04 2.92 



nosity function in IBlanton et al.l l|2003fl . 

From redshift layer 0.1 < z < 0.3 to 0.3 < z < 0.5, 
the estimated faint-end slope is found to c hange from 
a = -1.25 ± 0.03 to a = -1.35 ± 0.08. lEllis et alJ 
( 1996) constructed lminosity functions using a sample of 
1,700 galaxies observed by Autofib Redshift Survey out 
to z ~ 0.75 and gave that a steepens from —1.1 to —1.5 
with re dshift. Our resu lts is in good agreement with this 
trend. lLovedavl l|2004D studied the evolution of LF at 
z < 0.3 by a sample of 162,989 spectroscopic galaxies 
with magnitude limit of r < 17.6 in SDSS. The poorly 
constrained faint-end slope in redshift slice 0.2 < z < 0.3 
is due to the incompleteness in high redshift with a bright 
magnitude limit r < 17.6. Ioveday04 investigated the ef- 
fect of different absolute magnitude ranges on the esti- 
mated faint-end slope. It is found that the faint-end slope 
a changes from —1.17 to —2.18 if analysis of redshift slice 
0.1 < z < 0.15 is limited to Mo.i r < -21.5. This effect 
will propbably cause the overestimate of faint-end slope 
for redshift layer 0.3 < z < 0.5 in our sample: basically 
there are no low-luminposity points to tie down the faint 
end. Though not obvious of the incompleteness for the 
redshift slice 0.1 < z < 0.3 in Fig. 5, it could result from 
an incomplete sample, either. 

The estimated characteristic magnitudes M* = 
-21.80±0.16, M* = -21.91 ±0.12, and M* = -22.69± 
0.21 for the three redshift layers, are a bout 1 to 2 mag- 
nitude brighter than prev i ous results (lEllis et alJ 119961 
IBlanton et aIll200lL 120031 lLovedavl l2004[) . It is partly 
due to the different choice of the hubble constant with 
h = 0.75 with others h = 1. Another reason comes 
from the effect of larger magnitude errors and photomet- 
ric redshift errors than spectroscopic samples. We will 
demonstrate this effect in § 5. 

4. DEPENDENCE OF LUMINOSITY FUNCTION ON 
ENVIRONMENT 

To study the dependence of galaxy luminosity function 
on large-scale density environment, we subdivide galax- 
ies into high and low density two subsamples accord- 
ing to the density enhancements in three-dimensional 
redshift space. There are many methods implemented 
in the literature, here w e use the percolation algorithm 
l|Huchra fc GellerllrM^ . This algorithm can search iso- 
late group and cluster. The algorithm identifies every 
two galaxies by the projected separation D12 and the 
line-of-sight redshift separation zi2- 

D12 = sin(0/2)(zi + z 2 )c/H < D L {z u z 2 ,m u m 2 ), (8) 

z\2 = \zi - z 2 \ < z L {z 1 ,z 2l m 1 ,m 2 ) 1 (9) 
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creases from — 1.21±0.04 to — 1.37±0.04 and the charac- 
teristic magnitude brightens slightly from —22.26 ± 0.16 
to — 22.64±0.18. The right panel in Fig. 6 shows us the 
Fa (68.3% 2-parameter) , 2-cr (95% 2-parameter) and 3-a 
X 2 contours in the a — M* plane for high and low density 
populations. The difference is significantly obvious at the 
confidence level of 95%. Since the curves in Fig. 7 are 
normalized in the magnitude range —24 < Mj < —16, 
the bright-end density of underdense populations does 
not mean the numbers of luminous galaxies are larger 
than that in overdense regions. 



Fig. 6. — 1-cr, 2-cr and 3-a error contours for parameters M, and 
a. Left panel is that for the three redshift layers 0.03 < z < 0.1 
(solid line), 0.1 < z < 0.3 (dashed line) and 0.3 < z < 0.5(dot- 
dashed line) subsamples. Right panel is that for overdense (solid 
line) and underdense (dot-dashed line) subsamples. 



where z\ and zi refer to the redshifts of the two galaxies 
in the pair, 9 is their angular separation, and Dl and 
Zl are scaled to account for the magnitude limit of the 
galaxy catalog. All pairs linked by a common galaxy 
form a group. The limiting number density contrast is 



Sp 
P 



M, 



<$>{M)dM\ 



1. 



(10) 



where $(M) is luminosity function, Mi is the faintest 
absolute magnitude for galaxies with magnitude limit at 
fiducial distance. For a choosen density contrast and 
luminosity function (assumed without evolution for the 
simplicity), the critical distance Dq can be calculated. 
Taking into account the decrease of galaxy numbers with 
increasing distance, link parameter R can be calculated 
by: 



J^(M)dM 



R=[ 



J M ^ 2 $(M)dM' 



/3 



(11) 



In R, M\2 is the faintest absolute magnitude for a galaxy 
with magnitude limit at the mean distance of the two 
galaxies. Then Dl = RD$, and zl — Rzq. The magni- 
tude limit here is the complete magnitude. The fiducial 
redshift Zq we choose is the minimum redshift of the sam- 
ple, Zf = 0.03. Do is the fiducial distance corresponding 
to Zf. $(M) is chosen as that measured by Blanton et al. 
(2001) for local galaxies by SDSS commissioning data. 
The LF parameters are: <j>* — 1.46 x 10~ 2 /i 3 Mpc~ 3 , 
M* = —20.83, and a = —1.20. By this criteria, we can 
obtained our subsamples by Sp/p. Sp/p = 30 is chosen 
for our classification, with the fiducial distance D cor- 
responding to about 0.35 Mpc. In this scale, the high 
density subsample includes group galaxies and cluster 
galaxies (members larger than 5) and low density sub- 
sample contains field and void galaxies. We obtain two 
subsamples with 54,267 and 15,404 galaxies, respectively. 

Fig. 7 shows us the luminosity functions derived for 
subsamples classified by density environment. The dis- 
crepancy in fits may be due to the same reason given 
above. The best-fit Schechter parameters along with the 
number of galaxies considered in each density environ- 
ment are listed in Table 1. From Table 1 we can see that, 
with the enhancement of density, the faint-end slope in- 




6p/p>30 
Sp/p<30 



Fig. 7. — Luminosity functions fitted for overdense and under- 
dense subsamples. 



iCroton et alJ l)2004|) indicates the luminosity function 
by density environment. They implement galaxies ob- 
served by 2dFGRS with median survey depth is z ks 0.11. 
The local density contrast is determined by 5s- The same 
tendency of a is given for regions of different density 
contrast. In numerical simulations, Mo et al. 2004 stud- 
ied the dependence of the galaxy luminosity function on 
large-scale environment in hierarchical cosmology. The 
results predict that the characteristic luminosity, L* , in- 
creases moderately with density and the faint-end slope 
is quite independent of density in total, a is virtually 
constant for late types and increases from —1.3 in un- 
derdense regions to —1.8 in overdense regions for early 
types. Although our samples are not split into subsam- 
ples by galaxy types, the predictions are broadly in good 
agreement. The steepen of the faint-end slope from un- 
derdense environment to ove rdense environment is ex- 
plained in iTullv et all l|2002]) by a process of photoion- 
ization of the IGM which suppressed dwarf galaxy for- 
mation. Overdense regions typically collapse early and 
can form a dwarf galaxy before the epoch of reionization. 
Underdense regions, however, collapse later and are thus 
subject to the photoionization suppression of cooling 
baryons. The brighten of the characteristic magnitude in 
dense re gions is consistent with the morph ology-density 
relation Ipresslerl 119801 IBinggeli et al.l 11988^ . which is 
that the population in low density subsamples is dom- 
inated by late types and cluster regions have a relative 
excess of most luminous early-type galaxies. In structure 
evolution models, it is explained that, the most dense 
regions of the universe will have collapsed earlier, have 
larger merger rate and will contain more massive early 
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type galaxies. 

5. DISCUSSION AND CONCLUSION 

The uncertainties of galaxy photometric redshifts are 
about from 0.01 to 0.04 in different photometric magni- 
tude uncertainties. This is the best accuracy achieved by 
multi-color photometric information. This uncertainty, 
however, is much bigger than that of spectroscopic red- 
shift, which can be 0.0001. To calculate luminosity func- 
tion by photometric redshifts, we need to know how the 
rcdshift uncertainty and photometric magnitude uncer- 
tainty affect the shape of luminosity function. To as- 
sess this effect, we perform a series of simulations to 
fit luminosity function for a galaxy sample distributed 
as given luminosity function. 100 galaxy samples are 
created to evaluate the uncertainty of fitted M* and a. 
Each sample contains 5000 galaxies homogeneously dis- 
tributed in the 1 deg 2 cone-shaped comoving space in 
the redshift range z < 0.6. Gaussian distributed photo- 
metric errors with Am = 0.05, 0.20, 0.50 are assumed. 
Gaussian distribution with a Gaussian kernel of width 
(j z = a z (l+z) for photometric redshift erros are adopted, 
where a z is the redshift rm s residual at zero redshift 
1 I200T1 



<|Fernaridez-Soto et all 120011: iGhen et alJ I2003I) . o z = 
0.05, 0.20, and 0.50 are assumed. Here we choose a given 
luminosity function with M* = —21.17 and a — —1.26. 

Table 2 lists the estimated rms errors of M* and a 
for different photometric uncertainties and redshift un- 
certainties. We can find that, the faint-end slope a be- 
come flatter and the characteristic magnitude M* be- 
come brighter for all uncertainties. It means that the 
shapes of luminosity function become flatter, which is 
demonstrated intuitionisticly in Fig. 8. Fig. 8 shows 
us the luminosity function distribution in Gaussian dis- 
tributed photometric uncertainties of 0.05 and redshift 
uncertainties of 0.05 and 0.20. Fig. 9 gives the distribu- 
tion of Ms, and a in the 100 simulations with the same 
uncertainties as left panel in Fig. 9. From Table 2, with 
the increase of redshift uncertainties, the characteristic 
magnitude M*, is determined with an offset around 0.2 
mag brighter and an rms uncertainty about 0.1 mag, and 
a is determined with an offset around 0.2 flatter and an 
rms uncertainty about 0.02. The errors become larger 
with the increase of photometric and redshift uncertain- 
ties. We can see by simulation that photometric errors 
and redshift errors are two major factors for the mea- 
surement of luminosity function. This change trends of 
M* and a can explain partly the difference between our 
estimates and previous results. 

We study the evolution of galaxy luminosity function 
in the restfram r-band and the dependence of luminos- 
ity function on density environment out to z ~ 0.5 by 
69,671 galaxies composed by BATC sky survey and SDSS 
sky survey. We further the depth by photometric red- 
shifts and adopt the Schechter funtion as luminosity func- 
tion model. The evolution of galaxy luminosity func- 
tion is studied by three redshift layers 0.03 < z < 0.1, 
0.1 < z < 0.3 and 0.3 < z < 0.5. We subdivide the den- 
sity environment by criteria of isodensity contour. The 
density contrast is choosen to be Sp/p= 30. The principle 
conclusions are summarized as follows: 

(i) By simulation we find that photometric uncertainty 
and redshift uncertainty are two major factors that effect 
the measurement of luminosity function. In photomet- 




theoretlcal 
points with errors 




theoretical 
points with errors 



Fig. 8. — Luminosity function distribution of the 100 times 
simulations in the Gaussian distributed photometric uncertainty 
of 0.05 and redshift uncertainty of 0.05 and 0.20. 



M. = -21.17 
M. p = -21.37 
o„ = 0.076 





„.„- -1.26 
Q ,= 0.014 


r 


I " 



-1.3 -1.2 -1.1 -1 -0.S 



Fig. 9. — Distribution of fitted M* and a for the simulated 
samples same as that of left panel of Fig. 8. 



Table 2 

Simulated results of the accuracy of luminosity 
function parameters in different magnitude 
uncertaineties and redshift uncertainties. the 
given paramters of luminosity function are 
M, = -21.17 and a, = -1.26. 



Cm 




Af*. fitted 


°"M, 


fitted 




0.05 


0.05 


-21.28 


0.076 


-1.28 


0.014 


0.05 


0.20 


-21.53 


0.137 


-1.19 


0.017 


0.05 


0.50 


-21.38 


0.123 


-0.93 


0.022 


0.20 


0.05 


-21.28 


0.071 


-1.28 


0.013 


0.20 


0.20 


-21.53 


0.147 


-1.21 


0.028 


0.20 


0.50 


-21.38 


0.120 


-0.94 


0.023 


0.50 


0.05 


-21.28 


0.065 


-1.28 


0.015 


0.50 


0.20 


-21.53 


0.118 


-1.18 


0.017 


0.50 


0.50 


-21.48 


0.117 


-0.95 


0.028 



ric uncertainty of Am < 0.50 and redshift uncertainty of 
<j z < 0.50, the characteristic magnitude M», can be es- 
timated 0.2 mag brighter with a rms error about <Jm,= 
0.1; and the faint-end slope a, can be recovered with a 
factor 0.2 flatter and small uncertainty about a a = 0.02. 

(ii) There is slight evolution in the shape of galaxy lu- 
minosity function with observational depth. We further 
the accurate measurement of galaxy luminosity function 
to z < 0.3. The faint-end slope steepens slightly from 
— 1.21 to —1.25, —1.35 with the increase of redshift from 
0.03 < z < 0.1 to 0.1 < z < 0.3, 0.3 < z < 0.5. The 
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change of a is broadly c o nsistent wi t h pre vious claims 
such as lLlhV et all {1995^ : lEllis et all lfl99^ . 

(iii) Luminosity function differs distinctly with density 
environment. The faint-end slope for high density galax- 
ies is steeper than that for low density galaxies. The 
value of a changes from —1.21 at underdensity regions 
to —1.37 at high density regions. M* brightens from 
-22.26 to -22.64. 

We would like to thank the referee Dr. Jon Love- 
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